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Despite the growing literature on renewable energy sources, causal relationships between the variables 
that are selected as inputs of the models proposed in forecasting studies have not been investigated so 
far. In this paper, a novel approach to decide prediction input variables of wind and/or temperature 
forecasting models is suggested. This approach uses time series techniques; more specifically, Granger 
causality and impulse-response analyses between some meteorological variables. To conduct our study, 
wind speed, temperature and pressure data obtained from different regions of Turkey are employed. The 
results suggest that bidirectional causal relationships exist between these variables and that short-run 
dynamics differ with respect to location (inland versus coastal area). From this, it is concluded that 
renewable energy models must be built accordingly to improve prediction accuracy. 
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1. Introduction 


The usage of renewable sources such as solar and wind are 
exponentially increasing day by day. Accordingly, with the increasing 
attention being paid by energy industries, governments, non- 
governmental institutions or the general public to issues relating to 
the use and development of renewable energy technologies, the 


* Corresponding author at: University of Paris Ouest, Department of Economics, 
200 Avenue de la République, 92001 Nanterre cedex, France. 
Tel.: +33 140977815; fax: +33 140974198. 
E-mail addresses: hocaoglu@aku.edu.tr (F.O. Hocaoglu), 
fkaranfil@u-paris10.fr, fatihkaranfil@gmail.com, fatih_9@yahoo.com (F. Karanfil). 


1364-0321/$- see front matter © 2013 Elsevier Ltd. All rights reserved. 
http://dx.doi.org/10.1016/j.rser.2013.07.054 


research in this field has experienced a remarkable impetus in recent 
years. As will be discussed later, one of the avenues that has been 
investigated by researchers is the prediction of renewable energy 
generation. A major challenge that the scholars in this area have to 
face is that since meteorological variables carry stochastic behavior, it 
is not straightforward to determine which variables carry the most 
important information to predict another one. In other words, what- 
ever the prediction method employed the output of the prediction 
depends to some extent on some a priori assumptions on the 
relationship between the meteorological variables involved in the 
analysis. Consider the following example to clarify further this point. 
From a statistical viewpoint, in a clear region, temperature data and 
solar radiation data should be strongly correlated. However, in a 
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cloudy region, the correlation between the two is likely to be less 
apparent. In this perspective, it becomes clear that, for a specific 
region, if the existence of a possible relationship between meteorolo- 
gical variables is investigated and in the case of a relationship, if the 
nature of this relationship can be determined in terms of causality 
before a prediction model is developed and estimated, then both the 
accuracy of the model and the performance of the prediction can be 
improved. That is to say, a pre-analysis method should be developed 
and used before the prediction step. This is what the present paper 
aims to achieve. Therefore unlike the previous articles, in this study a 
pre-analysis methodology is suggested to test for the existence of a 
relationship between the methodological variables (both predicted 
variables and those used for prediction) and to explore the causal 
nature of this relationship. 

Essentially the method presented in this study consists of using 
time series techniques to assess the inter-relationship between a 
number of variables for renewable energy sources. In fact the use of 
such techniques has been receiving considerable attention in recent 
years especially with the increased focus on energy demand forecast- 
ing (see [1] for a recent review of this literature). The approach 
proposed in this paper differs from the exiting literature in that it 
focuses on both the possible causal chains among the meteorology- 
related variables and the magnitude of their deviations from the 
steady-state at each point in time. Our idea is that; a variable (for 
example, in our case, temperature or pressure) should be used for 
forecasting another (for example wind speed), only if there is a 
causality running from the former to the latter. In a seminal paper, 
Granger argues that this causality technique can be used for determin- 
ing whether one time series variable is useful in forecasting another 
and as discussed below, in its general form, it provides a cross- 
disciplinary framework for studying different relationships in different 
contexts [2]. 

It should be noted that the econometric methods used in this study 
are not new in the economics literature. However, it appears that the 
idea of using time series techniques with meteorological data is new 
in renewable energy literature. More specifically, the novelty of this 
paper is that it shows the way in which the information required for 
renewable prediction studies can be obtained using Granger causality 
procedure. The present analysis goes further by providing simulations 
of generalized impulse response functions (GIRFs) which makes it 
possible to trace out the transmission mechanisms revealed by the 
meteorological data used in the analysis. More precisely the GIRFs are 
used to determine the effect of a shock to any one of the variables on 
the current and future movements in others. Hence, this final analysis 
provides further evidence on the short-run dynamic linkages among 
these variables and how each of them responds to shocks by the other 
variables in the model. 

In view of the above, this paper contributes to two strands of the 
literature which have been unconnected so far. First, although its main 
objective is not to make a prediction, it still falls into the literature on 
renewable prediction since it provides an econometric procedure that 
enables a coherent model construction for predicting energy resource 
data and a better understanding of the dynamics of the relationships 
between the renewable energy sources. In other words, using the 
causality tests and GIRFs proposed herein, one would be able to decide 
ex ante on the variables that should be involved in the prediction 
analysis. Second, the present study adds to the growing body of 
empirical literature on Granger causality. Using time series analysis as 
a pre-analysis of renewable prediction, the paper enriches the 
application fields of Granger causality and GIRFs. 

The organization of the rest of the paper is as follows. Section 2 
gives a review of previous work on renewable energy prediction. 
Section 3 is divided into two parts: while in Section 3.1 the data used 
in the study are presented and discussed, in Section 3.2 the proposed 
method mentioned above is described and a very brief review of the 
econometric literature using time series techniques is included. The 


empirical results are illustrated and analyzed in detail in Section 4. 
Finally, the conclusion and future works are mentioned in Section 5. 


2. Literature review 


As mentioned in the previous section, accurate prediction 
models are of vital importance for energy related issues. Therefore 
a lot of prediction models developed, tested and published pre- 
viously. In this section deliberately selected ones are reviewed. 

Benghanem et al. developed an artificial neural network (ANN) 
for daily solar radiation modeling [3]. Soares et al. applied a 
perceptron neural-network (NN) technique to estimate hourly 
values of the diffuse solar radiation [4]. Azadeh et al. presented 
an integrated ANN approach for predicting solar global radiation 
by meteorological variables [5]. Cao and Lin proposed diagonal 
recurrent wavelet neural networks to forecast hourly and daily 
global solar irradiance [6]. In a similar study, Mellit et al. combined 
wavelet theory and neural networks and proposed a Wavelet- 
network model to predict daily total solar radiation [7]. In that 
study, various numbers of total solar radiation data were taken as 
inputs and different structures were developed. Hocaoglu et al. 
developed a novel hourly solar radiation forecasting methodology 
[8]. In a different study, a novel two dimensional modeling 
procedure is proposed [9]. Since photovoltaic sizing algorithms 
gives more accurate sizing results in case model generated data 
used as presented in [10], Hocaoglu developed novel analytical 
models for photovoltaic based system sizing algorithms and 
compared modeling performances with previously developed 
analytical models [11]. In another study, a different ANN based 
model was developed for the prediction of solar energy potential 
in Nigeria [12]. In that study meteorological geographical inputs 
were used in prediction. Behrang et al. used daily mean air 
temperature, relative humidity, sunshine hours, evaporation, and 
wind speed values to predict daily global solar radiations [13]. 
According to an extensive and in-depth review of the literature on 
the use of artificial intelligence (AI) in photovoltaic system 
modeling by Mellit et al., Al techniques provide the possibility 
for sizing PV-systems with reasonable accuracy even in the case of 
a lack of complete data [14]. A more general review of solar energy 
modeling and prediction techniques can be found in Khatib et al. 
[15] and also in Besharat et al. [16]. 

On the other hand, for wind speed prediction, Salcedo-Sanz 
et al. proposed a method to improve the accuracy of the wind 
speed prediction systems based on exploiting diversity in the 
input data of the neural networks [17]. Kavasseri and Seethara- 
man used fractional f-Autoregressive Integrated Moving Average 
(f-ARIMA) models to predict day-ahead wind speed forecast [18]. 
Louka et al. used Kalman filtering to improve wind speed forecast 
accuracy [19]. In another study, a history based method is devel- 
oped to predict short-term ahead wind speed [20]. Mohandes 
et al. used support vector machines in the same topic [21]. Grassi 
and Vecchio suggested a two-hidden layer neural network to 
predict the wind energy output [22]. From the same perspective 
Haque et al. proposed soft computing models (SCMs) “augmented” 
by a similar days (SD) method which is found to lead to an increase 
in the level of performance of wind speed forecasting technique 
[23] (see also [24] for a review of the use of ANN in this field). In a 
different study, temperature values are employed as model inputs 
and accurate wind prediction results are obtained using hidden 
Markov models [25]. Cadenas and Rivera compared ARIMA and 
ANN models in wind speed predictions [26]. 

In most of the previous studies, the prediction model inputs are 
decided intuitively. In this work, the necessity of a pre-analysis for 
prediction models is emphasized and a novel approach is sug- 
gested. This approach consists of using a time series analysis on 
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Fig. 1. Temperature, wind speed and air pressure data measured at Antalya region. 


some meteorological variables in order to investigate time series 
properties of each of these variables and then to search for a causal 
relationship among them. This approach can provide clear indica- 
tion of which variables should be incorporated in any prediction 
model as an input, and in consequence, more accurate prediction 
models can be developed. 


3. Data and methodological issues 
3.1. Data description 


To develop the econometric model, which will be discussed in the 
following section, hourly measured yearly wind speed, temperature 
and air pressure data obtained from Antalya, Izmir and Kayseri 
regions are selected. To prevent the developed approach from 
possible errors, the measured data are first pre-analyzed. This pre- 
analysis includes the followings: the missing data that are absent in 
the whole set are interpolated using future and previous data. Hence, 
erroneous data with invalid values were rectified as a pre-analysis. 

The wind speed, temperature and air pressure data are illu- 
strated in Figs. 1-3 for Antalya, Izmir and Kayseri regions, 
respectively. The range of data is shown in Table 1. 

It is obvious from Figs. 1-3 and Table 1 that, not only the ranges 
but also the variations of the data are different. Hence, the 
deliberately selected data has different characterizations. 


3.2. Methodology 


In recent years, the use of Granger causality procedure have 
gained increased popularity and today it is performed frequently 
in a widespread field of applications such as economics or epilepsy 
analysis.' In some promising studies, Granger causality tests are 
employed to increase our understanding of climate variability. 
To give some examples of the previous works close to our own, 
Mosedale et al. provide an application of Granger causality in an 
ocean-atmosphere system and study the influence of daily sea 
surface temperatures (SSTs) on daily values of the North Atlantic 
Oscillation (NAO) [30]. They conclude that there exists a significant 
feedback of SSTs on the NAO. In another interesting and related 
paper, Smirnov and Mokhov investigate the effects of variations in 
carbon dioxide atmospheric content, solar activity, and volcanic 


1 Here are some examples for the use of the Granger method in different areas 
of research. The most related literature to the present study using Granger causality 
tests investigates causal relationships between environmental variables (e.g. CO 
emissions) or energy variables (e.g. primary energy consumption) and economic 
development indicators (such as gross domestic product (GDP) or industrial value 
added). See Ref. [27] for an overview of the literature in this field. On the other 
hand, in the more distant literature, in epilepsy studies for instance, Granger 
causality tests are generally used for determining direction of relationships 
between local field potentials (see for an example [28]) or for examining structures 
of different neurological states (see for example Ref. [29]). 


EO. Hocaoglu, F. Karanfil / Renewable and Sustainable Energy Reviews 28 (2013) 204-214 207 


p (mb) 
1030 4 
1020 - 
1010 4 
1000 - 
990 + 
ss 
0 2000 4000 6000 8000 
hour 
t (Celsius) 
40 


30 


20 


10 


0 2000 


w (m/s) 
15 
10 
5 
0 
T T T T T 
0 2000 4000 6000 8000 
hour 


4000 6000 8000 
hour 


Fig. 2. Temperature, wind speed and air pressure data measured at Izmir region. 


activity on the global surface temperature (GST) [31]. They show 
that all the three factors Granger cause the GST. 

Since we have given a description of the variables involved and 
the basic structure of our idea, we now proceed with our econo- 
metric analysis which consists of testing for causal relationships 
and dynamics between the variables of interest, namely air 
pressure (p), temperature (t) and wind speed (w). It is well known 
that the order of integration of all variables used in the analysis 
has to be established before deciding on which time series method 
to use. For example, if the series are found to be stationary, the 
standard Granger causality test in an unrestricted vector autore- 
gressive (VAR) model can be conducted by testing whether 
parameters of explanatory variables are jointly zero. On the other 
hand, in the same framework, non-stationarity series having a unit 
root is a random walk and cannot be used for forecast. This 
problem has given rise to alternative testing procedures such as 
error-correction model (ECM) [32,33] and autoregressive distrib- 
uted lags (ARDL) bounds test [34]. 

In line with the existing literature on time series analysis, first 
of all, time series properties of the data are checked by performing 


? Detailed presentation of the time series analysis and Granger causality 
procedure can be found in Hamilton [35, Chapters 11 and 19]. For a brief discussion 
on the same issue, see for instance [36]. 


the augmented Dickey-Fuller (ADF; [37]) and the Phillips and 
Perron (PP; [38]) unit root tests based on the model given by the 
following equation: 


k 
AX; = ao + But 01Xt-1 + È AAXi + Ur ad) 
i=1 


where X is the variable to be tested, y is the trend variable, A is the 
first-difference operator and us is Gaussian white noise. On the 
other hand, £ is the coefficient on the time trend and k is the lag 
order of the autoregressive process. For any given time series, if 
a, = 0 then this series can be characterized as unit root or random 
walk. In other words, in this case the variable behaves as a non- 
stationary process. As mentioned above, in order to address the 
non-stationarity issue, several methods have been suggested in 
the literature. However, if all variables are found to be stationary, 
this result will bring us the standard Granger causality framework, 
since it is not necessary to use more complicated models 
Granger gives the definition of causality in terms of prediction: 
given to time series x and y, “if y, causes x;, then x;,; is better 
forecast if the information in y,_; is used than if it is not used” 
[39, p. 200]. In time series analysis, this means that one may 
determine whether one time series (here the variable y) is useful 
in forecasting another (here the variable x). In our case, this 
statement implies that the Granger causality procedure can be 
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Fig. 3. Temperature, wind speed and air pressure data measured at Kayseri region. 


Table 1 
Range of the data. 


Regions Temperature (°C) Wind speed (m/s) Air pressure (mb) 
Antalya 1.3-40.06 0-15.3 994-1020 
Izmir —0.8-38.2 0-13.6 995-1028 
Kayseri —16.4-37.5 0-11.5 873-903 
Table 2 


Granger causality test hypotheses. 


Direction of causality Test thypotheses 


Xı does not cause X2 Ho: bj; =0 Vi=1,...,m 


Hı : 4b4,;40 Vi=1,...,m 


Xı does not cause X2 Ho: @;=0 Vi=1,...,n 
Hi: Jaz 0 Vi=1,...,n 


well suited for wind speed, temperature and pressure forecasting. 
Before we go on to discuss the econometric methodology, to avoid 
any possible confusion due to the use of the term “prediction” 
herein, we would like to underline the following point, although it 


Table 3 
VAR Granger causality test results (Block exogeneity Wald test). 


Independent variable: 
air pressure (p) 


Independent variable: 
wind speed (w) 


Independent variable: 
temperature (t) 
Direction of Chi- 


Direction of Chi- Direction of Chi- 


causality square causality square causality square 
Estimation results for Antalya region (included lag number: 4) 

wp 65.3* pow 60.2* pt 230.6* 
top 366.9% tow 175.9% w-t 99.5* 

Estimation results for Izmir region (included lag number: 7) 

w>p 58.9* pow 89.3* pot 256.5* 
top 518.7* tow 479.7" wat 192.3* 
Estimation results for Kayseri region (included lag number: 6) 

wp 38.2* pow 122.8% pot 100.8* 
top 564.6" tow 297.3* w-—t 123.6° 


* Significance at 1% level. 


has been already discussed above and will be emphasized again in 
the following sections: Granger's method for detecting causality 
between two or more variables cannot be used to predict any 
variable. Instead, indicating the direction of causality (if some 
causal relationships exist), it makes it possible to ascertain which 
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variable may be used to predict the other. In consequence, the 
procedure proposed in this paper should be considered as a pre- 
analysis method for renewable prediction. The variables among 
which causal relationships are detected may be used ex post to 
predict energy resource data. 

Technically speaking one can show that the process X; does not 
Granger cause the process X, if 


E(X2t/lt-1 (X2), l-1 (X1) = E(Xae/Tt-1(X2)) (2) 


where /;_1(X;) is the space generated by the linear combinations of 
the past values of Xj. 

With stationary series, following Granger, this can be done by 
performing a set of tests based on the regressions given below [2]: 


m n 
Xy=at DY auXyit DY aaiXo-i+ ete 
i=1 i=1 


m n (3) 
Xx =h+ Y biiXii + Y b2iXat-i + ert 
(24 


isi 


where a and £ are constant terms, ¢;; and epp are white-noise 
series and m and n represent the lag orders. The classical Fisher 
test can be used to detect the causal chains between X, and X3. 
Granger causality test is then formulated as the hypothesis 
problems given in Table 2. 

We should also mention that this bivariate framework can be 
extended including more variables in the regression equations. 
Moreover, it is better to do so, since, as Lutkepohl argued in [40], 


Granger causality test in a bivariate framework may be subject to 
the omitted variables bias. Multivariate systems become widely 
used in some recent studies on time series econometrics. In our 
model we use temperature, wind speed and, air pressure in order 
to detect all possible causal chains between them. The results are 
given in the following section. 


4. Empirical results 


As we have discussed above, the first step of our empirical 
analysis consists of testing for the stationarity of the variables 
involved in our study. Using Eq. (1), ADF and PP unit root tests 
suggest that all series are stationary, in other words, they are all 
integrated of order 0, that is I(0) (see Appendix A for the results of 
unit root tests). In this case, the standard Granger causality test 
should be employed in a multivariate VAR framework and the 
system of equations given in Eq. (3) evolves to the following 
system: 


a b c 
P= + 2 PriiPt-i + x By2iWe-i + x Aisitti + Ure (4.1) 
t ng t ms i= 
a b È 
We=W2t+ È PaiPit È Boaiweit È Positri + uae (4.2) 
iSi iZi izi 


Response to Generalized One S.D.Innovations + 2S.E. 


Response of AN_P to AN_P 


Response of AN_P to AN_T 


Response of AN_P to AN_W 


0.6 -~—___ 0.6 0.6 
0.2 + 0.2 4 0.2 ees 
0.0 0.0 He Ys aes ae 
024 of Seen oe | 
-0.4 : i : : -0.4 , f i i -0.4 : ; : i 
5 10 15 20 5 10 15 20 5 10 15 20 
ae Response of AN_T to AN_P oe Response of AN_T to AN_W 
2.0 4 2.0 -| 
1.5 - 15 J 
1.0 - 1.0 4 
0.5 - 
0.0 paeee == 
l | | l e ee 
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Response of AN_W to AN_P 
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Fig. 4. Generalized impulse-response functions for Antalya. 
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Response to Generalized One S.D. Innovations + 2S.E. 
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Fig. 5. Generalized impulse-response functions for Izmir. 


a b € 
te=y3 + $ P3iPtit È P32Wi-i + © P33itt-i + Use (4.3) 
i=1 = iZi 


where all notations (depending on where they are used) are 
similar to those in Eq. (3). Lag orders a, b and c are selected based 
on the lag exclusion tests. Significance of the parameters f; for 
each variable in each equation will indicate the direction of 
causality. 

We can now apply the standard Wald test for this purpose. The 
results are reported in Table 3. 

From Table 3, it can be seen that: (1) all parameters are found 
to be significant at the 1% significance level; (2) the neutrality 
hypothesis (no causal relationship exists) can be rejected at the 
1% level; (3) this result holds for all three regions studied and all 
variables involved without any exception. More specifically, the 
significance of all parameters in the system of equations (given by 
Eqs. (4.1)-(4.3)) implies that there exists bidirectional causality 
among these variables.’ Taking all together, these findings indicate 
that if one attempts to conduct a renewable prediction study, 
temperature, air pressure and wind speed can all be used to do so. 

The final step of our empirical work consists of plotting GIRFs 
of each of the variables for each of the cities. In fact, this method, 
first introduced by Koop et al. [41] and Pesaran and Shin [42], 
enables one to assess how each variable responds to innovations in 
other variables, and to determine whether the shock is permanent 


3 We have conducted the same analysis including solar radiation instead of 
temperature and found identical results to those reported here (see Appendix B). 


or not. Figs. 4-6 show the responses of pressure, temperature, 
and wind speed to one standard deviation shocks to other 
variables in the VAR system, including the plus and minus two 
standard deviations from the mean (representing thus the 95% 
confidence level). 

Figs. 4-6 present the impulse-response paths of the variables 
involved in the VAR model up to 24h after a one standard 
deviation shock is stimulated from others. Overall, comparing 
the GIRFs for each of the three cities we find very similar results 
with some differences in sensitivity. It follows also from these 
figures that there is a considerable self-response (i.e. response to 
own shock) for all the three variables in the model. Only for the 
case of wind speed, although it increases in response to its own 
shock, it moves smoothly downward to its pre-shock level. The 
oscillating trajectories between different variables when tempera- 
ture is considered come from the fact that the GIRFs are computed 
over the entire 24-hour period. However, note also that the level 
and the period of these oscillations are different with respect to 
impulse-response variables. For instance, in Antalya (Kayseri) 
temperature tends to have positive (negative) responses to the 
shock of pressure at the initial period and then oscillates about 
zero (dies out). On the other hand, during the 24-hour period, the 
response of wind speed to a positive one standard deviation shock 


4 We prefer to use GIRFs since this method has the advantage over traditional 
impulse response analyses in that it is invariant to the ordering of the variables in 
the VAR model. 
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Fig. 6. Generalized impulse-response functions for Kayseri. 


Table A1 
Results of unit root tests for air pressure, temperature and wind speed. 


Table B2 
Causality test results for the VAR system involving radiation instead of 
temperature. 


Independent variable: 


Variables Augmented Dickey-Fuller (ADF) test Phillips-Perron (PP) test ; : 
es ee Independent variable: Independent variable: 
Antalya Izmir Kayseri Antalya Izmir Kayseri air pressure (p) wind speed (w) 

p -5.47 -540 _753 -5.04 -5.66 —681 Direction of Chi- Direction of Chi- 

t -239 -3.27 -3.02 -6.39 -396 -6.43 causality square causality square 

w — 11.48 — 12.14 — 12.83 —22.76 —25.53 —26.48 


Critical values for the 1%, 5% and 10% levels of significance are —3.43, — 2.86 and 
—2.56, respectively. Lag determination for the ADF test is based on the Akaike 
Information Criterion (AIC). 


Table B1 
Results of unit root tests for radiation. 


Variable Augmented Dickey-Fuller (ADF) test Phillips-Perron (PP) test * Significance at 1% level. 
Antalya Izmir Kayseri Antalya Izmir Kayseri 
r —5.45 — 5.02 — 5.76 —13.79 —13.72 -14.78 


Critical values for the 1%, 5% and 10% levels of significance are — 3.43, — 2.86 and 
—2.56, respectively. Lag determination for the ADF test is based on the Akaike 
Information Criterion (AIC). 


radiation (r) 


Direction of 
causality 


Estimation results for Antalya region (included lag number: 7) 


w>p 
r>p 


107.1* 
654.4* 


p>w 
row 
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por 
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Estimation results for Izmir region (included lag number: 5) 


por 
w-r 


por 


wp 23.1" pow 37,2" 

r>p 827.4* row 534.4* 

Estimation results for Kayseri region (included lag number: 6) 
w>p 62.3* pow 174.5* 

r>p 985.8* r>w 242.2* 


w>r 


Chi- 
square 


937.9" 
62.4* 


435.1* 
724" 


338.1* 
63.0* 


in the pressure level never becomes positive for the case of 
Kayseri, contrary to the cases of Antalya and Izmir, which are 
coastal regions. Another difference that can be observed between 
coastal and inland regions is that wind speed tends to have a 
positive impact on the temperature in Antalya and Izmir while the 
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most evident negative impact of wind speed on temperature is 
found to be in Kayseri, which is an inland region. These findings 
indicate that regional differences may lead to different short-run 
dynamics and impulse-response characteristics. Finally it should 
be added that the estimated responses to all shocks are transitory 
and take 40-240 h to die out (see Fig. C1 in Appendix C) and that 
of the nine shocks for each of the cities, the self-response of 
pressure is estimated to be the most persistent. 

The previous literature on renewable energy prediction has not 
examined the causal chains between the variables before con- 
ducting a prediction study. Among the studies examined in this 
paper, for instance in [5] seven different variables were considered 
as inputs of the model to predict solar global radiations. In [13] on 
the other hand, six different combinations of the input variables 
were employed to predict global solar radiations. However in 
these and similar studies, prior to model construction no causality 
analysis has been undertaken. It is evident that such a task 
proposed in this study will improve the accuracy of the forecast. 


Response of AN_P to Generalized One 


Response of IZ_P to Generalized One 
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Furthermore it may be possible to assign some weights to the 
variables of the model according to their impulse-response func- 
tions in order to develop advanced prediction models. Such 
analyses should be considered as future tasks of research in 
this field. 


5. Conclusions and future work 


In this study, a novel methodology is developed for renewable 
energy analysis. This methodology suggests a pre-analysis using 
Granger causality and GIRFs before construction of a prediction 
model. Using this methodology, the inputs of the prediction model 
should be determined and unnecessary inputs should be removed 
in order not to affect negatively the accuracy of the prediction 
model. In other words, the prediction performance and accuracy 
can be improved by taking into account the direction of causality 
among the variables that are to be introduced into the prediction 


Response of KA_P to Generalized One 
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Fig. C1. Generalized impulse-responses up to 240 h. 
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model. According to our results, bidirectional causal relationships 
can be established. Hence, taking into account differences in 
the impulse-response functions, we conclude that although it is 
possible to use all meteorological variables, in at least these 
three cities, regional differences (coastal versus inland region) 
may influence short-run dynamics between the meteorological 
variables which should be, thus, considered in any prediction 
study. This confirms our initial suggestion that the accurate 
variables that should be incorporated in the model must be 
selected using the suggested methodology. 

Future research should take into account regional disparities or 
differences in periods and intervals of time. Recall that this paper 
uses hourly data for one year, whereas it is possible that seasonal 
effects play a strong role in the relationship between the meteor- 
ological variables. Accordingly, in different seasons, different 
causal linkages and short-run dynamics may be found, implying 
that in prediction studies, different sets of variables may be used 
in different time intervals. We have not addressed the issue of 
seasonality in this paper which may prove fruitful as a future 
research topic. 
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Appendix A. Unit root tests 


Table A1 presents the ADF and PP test statistics for the levels of 
the variables with the null hypothesis of non-stationarity. 

The ADF test results indicate that all the series involved in 
the analysis are stationary, with the exception of the temperature 
variable t in Antalya. However, the PP test shows that all variables, 
including t in Antalya, are stationary even at the 1% level. In fact, it 
is commonly suggested that using the PP unit root test has several 
advantages over the ADF test since it does not require a lag length 
selection procedure and that the error term in the test equation is 
corrected for both serial correlation and heteroscedasticity [38]. 


Appendix B. Estimates of the VAR model including radiation 
series instead of temperature series 


The purpose of this appendix is to show that when we 
substitute temperature series (t) by solar radiation series (r) in 
the initially estimated multivariate VAR model (see Eqs. (4.1 )-(4.3)), 
the causality tests do not lead to different conclusions. For this, 
following the procedure already explained in the methodology 
section, first of all, unit root tests should be conducted to determine 
the order of integration of the variable r (Table B1). 

As was the result for the previous unit root tests for the 
variables w, p and t, both the ADF and PP tests indicate that the 
variable r is also a stationary series. Then the causality tests can be 
carried out by a VAR model composed of three variables (air 
pressure (p), wind speed (w) and solar radiation (r) instead of 
temperature (t)). The results are reported in Table B2. 

Form Table B2 it follows that bidirectional causal linkages exist 
between these three variables, as it is the case for the air pressure- 
wind speed-temperature nexus. 


Appendix C. The GIRFs of the regions up to 240h 


See Fig. C1. 
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